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THEORETICAL SPECTRA OF THE AM CVN BINARY 
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Taking into account a range of parameters determined from the evolution- 
ary models and available observational data, the detailed non-LTE spectrum 



p |. for the primary star and the irradiated donor star in the AM CVn system SDSS 

Q . J0926+3624 are constructed based on the TLUSTY stellar atmosphere code. The 

-)— > , combined spectrum of the primary and the donor stars along with a multi-color 

ci '. blackbody spectrum of the accretion disk that reproduces a detailed numerical 

model is compared to the SDSS optical spectrum of the system. The photometric 
Kj^ , flux of the primary star inferred from eclipse observations is compared with the 

"^ ! synthetic spectrum. The model fit of the two independent observations provides 

Ti^j- ! an upper limit on the distance of the system for different effective temperatures of 

the primary. In addition, an upper limit on the combined flux of the disk and the 
donor in the infrared region wherein the contribution of the primary is negligible 
is also determined. It is shown that the spectrum of a sufficiently cool donor can 
exhibit emission lines due to irradiation from a hot primary and the emission 
features should be detectable in the infrared even though the contribution of the 
flux from the disk dominates. Thus, it is pointed out that infrared observations 
5^ I of the system would provide important information on the thermal state of the 

donor as well as useful insight on the thermal properties of the primary star and 
the accretion disk. 
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1. INTRODUCTION 



The only known ec lipsing AM Canum Venaticorum object to date, SDSS J0926+3624 
(jAnderson et al. II2005I ). belongs to the class of ultra compact binary systems character- 
ized by a pair of interacting white dwarfs with the shortest orbital period among any 
other binary subclass. Such a system is driven by the loss of orbital angular momen- 
tum by the emission of gravitational waves and is among a class of systems thought to 
constitute the largest fraction of interacting ultra compact binaries in the Galaxy [see 
Nelemans. Yungelson Sz Portgies Zwart I fJ200ll . |2004| )]. An understanding of their properties 
is important for their detectability not only as individua l candidate sources, but also as back- 
groun d sources of low frequency gravitational waves [e.g. jNelemans. Yungelson fc Portgies Zwart 
J2OO4J ): Istroeer. Vecchio fc Nelemans I feoosh ]. 



A detailed understanding of the formation of SDSS J0926+3624 (hereafter J0926+3624) 
also allows one to place important constraints on the efficiency of the common envelope phase 
of binary evolutio n in forming short period binary systems containing double white dwarf 
components [e.g., iNelemans fc Tout I ( 120051 ): for a review of the common envelope phase of 



evolut ion for the general binary population, see llben fc Livio I (Il993[ ) ; iTaam fc Sandquist 



( 120001 )]. For AM CV n systems, in general, their formation may in volve a white dwarf or 
helium star channel ( iNelemans. Yungelson &: Portgies Zwart 1 120011). or an evolved hydro- 



gen main seq u ence star channel (jPodsiadlowski. Han &: Rappaport II2003I ). As shown by 
Deloye et al. I ( 120071 ) their subsequent evolution is affected by the thermal state of the donor 



star. 

The observational investigation of these double white dwarf systems provides important 
insights into their nature and formation. For example, optical observations of the line spec- 
trum from the system carry the combined physical and chemical properties of the accretor 
and the accretion disk while infrared observations of the spectrum may provide important 
constraints on the thermal state, chemical composition, and surface gravity of the donor. 
As a result, atmospheric studies of the donor in combination with knowledge of its distance 
can be used as an additiona l probe of the formation channel and evolutionary state of the 
system (JRoelofs et al. 1120071 ). 
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J0926+3624 offers a unique opportunity for probing the thermal properties of the donor 
star since its surface gravity, essential for calculating the synthetic spectrum, is known. In 
addition, the orbital separation between the accretor and the donor, which determines the 
amount of irradiation onto the donor, is known from the component masses and orbital 
period. 

In this Letter, we present the model composite spectra consisting of the non-LTE spectra 
of an irradiated donor, and a primary star as well as a multi-color blackbody spectrum of 
the accretion disk. The model parameters for the system are discussed in §2. The formalism 
for estimating the contribution of the accretion disk is also discussed in this section. In §3, a 
brief description of the numerical methods employed for calculating the non-LTE theoretical 
spectra of both the accretor and the donor stars is provided. The results are described and 
discussed in §4 followed by our conclusions in the last section. 



ATMOSPHERIC MODEL PARAMETERS 



The range for most of the model parameters of J0926+3 624 is constrained by t he mod- 
elling of eclipse light cu rves. Based on the theoretic al work of iBildsten et al. I ( l2006l ) and the 



observational study of iCopperwheat et al. I (120111 ) the mass range of the primary (donor) 
is estimated as 0.79 — 0.89 Mq (0.027 — 0.038 Mq). Taking the mass of the primary and 
the donor stars as 0.84 M^t^ an d 0.029 M© respectively and the observed orbital period of 



28.3 min, iDeloye et al. I (120071 ) estimate the radius of the primary to be 0.01 R© and the 
radius of the donor to be 0.043 Rq. Hence, we set the surface gravity of the primary at 
\og{g) = 8.3 and that of the donor at log{g) = 5.6 where the surface gravity g is in cm s~^. 
From the center-to-center distance and the radii of the components we find the surface-to- 
surface separation to be 1.67 x 10^° cm. For irradiation from the photosphere of the primary 
star onto the surface of the donor star, the surface-to-surface separation is considered in 
our calculations. Although the incident radiation impinges on one face of the donor, we 
consider complete redistribution of the irradiated heat energy over the entire surface of the 
donor. As the donor star is likely to be tidally locked with the primary star, its spin rotation 
period should be the same as the orbital period. The small orbital period and the relative 
sizes of the components would tend to establish a complete redistribution of irradiated en- 
ergy over the entire atmosphere of the donor star. The flux difference of the day and night 
sides of Jupiter size exoplanets or biting with a period of a few days ranges from 10% to 50% 
( iBurrows. Ibgui fc Hubeny II2008I ). In the case of an irradiated but self-luminous object with 
radius about an order of magnitude less than Jupiter and orbiting with a period as small as 
28 mins, the difference in the irradiated and non-irradiated sides should be much less. 
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The effective temperature of botli tlie primary and tlie donor remains u ncertain owing to 
the ab sence of a direct measurement of the distance to the system. Recently, ICopperwheat et al. 
( 120111 ) estimated the contribution of the primary white dwarf star from the echpse depth, 
measured from multi-band, high time resolution light curves of the system. By fitting the flux 
in the ultraviolet and optical bands with a synthetic spectrum they establish a lower limit to 



the ef fectiv e temperature o f the p rimary Tgff^p as 17000K. In previous work, iDeloye et al. 



(120071 ) and iBildsten et al. I ( l2006l ) theoretically estimated T^ffp to lie between 21000K and 
24000K based on its cooling and compressional luminosity. iDeloye et al. I (120071 ) also predict 
the effective temperature of the donor stai Tpjf^r, to lie in the range between 1750K and 
4500K. Guided by the blackbody model of IDeloye et al. I (120071 ). we adopt three different 
%ff^P (18000K, 21000K and 25000K) and three different %ff^D (3800K, 4000K and 4400K). 



The evolutionary model calculations by IDeloye et al. I (120071 ) estimate the accretion 
rate to be about 10~^° y\.p)'iir~^. There have been a few attempts to model the flux from 



the accretion disk (see lNagel. Ranch fc Werner I (120091 ) and references therein). However, all 
these model s assume a steady st a te dis k. Further, from the variation in the location of a 
bright spot, ICopperwheat et al. I (120111 ) found that the geometry of the outer disk departs 
from circular motion. A detailed non-LTE model of a non-steady state disk with non-circular 
geometry is beyond the scope of the present investigation. Therefore, we approximate the 
flux of the disk by considering 



disk 



Hou 



2nrB{r)dr 



Ri; 



where B{r) is the Pl anck function corresponding to the radial distribution of the effective 
temperature given by iNagel. Ranch fc Werner I (120091 ) 



T(r) 



SGMpM 




1/4 



(2) 



where Mp and Rp are the mass and the radius of the primar y respectively, M is the rnass ac- 
cretion rate and a is Stefan-Boltzmann constant. Following iNagel. Ranch fc Werner I (|2009[ ) 
we take the inner radius Rin of the disk equal to 1.3Rp. The outer radius Rout is taken to be 
0.45a which is the maximum value of the radius of the disk determined from the observation 



of lCopperwheat et al. I (12011! ) where a is the orbital separation. A c omposite blackbody spec 



t rum reproduces the non-LTE steady-state accretion disk model of iNagel. Ranch fc Werner 



(120091 ) well. The contribution of the disk to the total flux depe nds on the inclination angl e 
and for the present object the inclinati on angle is as high as 83° (ICopperwheat et al. Il201ll ). 
A comparison of the model spectra of iNagel. Ranch Sz Werner I (120091 ) for different inclina- 
tion angles implies that the flux should be reduced by a factor of 4 to 7 when the inclination 
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angle is 83°. Therefore, in order to obtain an upper limit on the contribution of the disk to 
the total flux, we assume that at 83° inclination angle the flux is reduced by a factor of 4 as 
compared to its value at a face on view. 



NON-LTE MODEL ATMOSPHERES 



The stellar-atmosphere code TLUSTY developed by Hubeny and colleagu es (IHubenv 



1988 



Hubeny fc Lanz Ill995l ) and the related spectrum synthesis code SYNSPEC ( Hubeny. Hummer fc Lan 



19941 ) are used for the investigation. In using this code, we have adopted the atomi c data 
and model atoms provided on the TLUSTY homepage! (JLanz fc Hubenv 1120071. 120031). T he 
atomic spectral line database from CD ROM 23 of R. L. Kurucz (JKurucz &: Bell Ill995l ) is 
used in SYNSPEC. The elemental abundances adopted are listed in Table 1. The elements 
whose abundance is less than 10~^ are ignored as they do not significantly contribute to the 
total opacity. We have incl uded micro-turbulence at the level of 2 — 5 kms~^ relevant for 
high effective temperatures (JLanz fc Hubeny 1120071 ). However, the overall spectrum is not 
altered even if this value is increased by a few times. For all the cases, the emergent flux 
is integrated over all frequencies and the energy balance was checked for self consistency 
corresponding to the effective temperature for models with and without irradiation. 

In order to achieve convergent solutions for the atmospheric model of the donor star, 
we have neglected the effects of convectio n. However, comparison with a model calcu lated 
by using the atmosphere code PHOENIX Barman. Hauschildt. Short fc Baron I (J2000l ) and 
references therein] implies that convection affects the spectrum in the optical where the flux 
of the donor is negligible as compared to that of the primary or the disk. 

TLUSTY is well te sted and widely used for constructing the sp ectrum of both hot stars 

( Lanz fc Hubeny 1120071 ) and ultra-cool objects such as brown dwarfs fJBurrows. Sudarsky fc Hubeny 
20061 ). We compared our results with a set of similar models, but with different elemental 
abundances constructed by using the atmosphere code PHOENIX. These models were kindly 
provided to us by F. Allard and D. Homier (private communication 2008). 



4. RESULTS 

Since the abundances of elements other than He are low, the opacity of the atmosphere of 
the donor as well as the primary is mainly determined by this constituent. Fig. 1 illustrates 



^ http://nova.astro.uind.edu/indcx. html 
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the thermal structure of the atmosphere of the donor star as a function of pressure for 
different levels of irradiation from the accretor. It can be seen that at a pressure of about 
100 bar, the temperature increases rapidly. This rapid temperature increase occurs as the 
Rosseland mean optical depth exceeds unity. The total pressure varies slowly with respect 
to the optical depth as the optical depth increases from 1.2 at about 576 bar to 116 at 
about 667 bar for T^ff^D = 3800K. When irradiation is included, the temperature varies 
significantly at optical depths less than unity. For an effective temperature as low as 18000K 
for the primary and as high as 4400K for the donor, a few emission lines appear only in 
the infrared as can be seen from Fig. 2 (panel A). As the irradiation energy increases with 
the increase in the effective temperature of the accretor, the spectrum changes dramatically 
with the appearance of emission lines in both the optical and infrared wavelength regions 
(see panel B in Fig. 2). For a fixed orbital separation, the emission features in the spectrum 
of the donor not only depend on the effective temperature of the accretor, but also on that 
of the donor itself. As can be seen in panel B of Fig. 2, weak emission lines are apparent 
from the optical, e.g., at wavelength slightly greater than 0.5 /im to infrared if the donor is 
cooler {Tcff^D = 3800K) whereas strong emission lines are present in the entire spectrum if 
Teff^P is increased to 21000K and to 25000K even if the donor is hotter than 3800K. 

The combined spectrum of the system is derived by multiplying the surface area of each 
component with its flux and then summing them up with the flux of the disk. In Figs. 3 and 
4 we present the synthetic spectra of the primary and the irradiated donor, the composite 
blackbody spectrum of the accretion disk as well as the combined flux which is compared with 
the well calibrated and phase-averaged Sloan optical spectrum of the object. The observed 
spectrum from 0.38 to 0.92 /xm is reproduced from the SDSS archival data. The flux of the 
primary at u', g' and r' bands are also presented for a comparison. We have not performed 
any statistical analysis while comparing the synthetic spectra with the observed flux and our 
fit is by eye. In Fig. 3, we also compare our synthetic spectra with that constructed by using 
the atmosphere code PHOENIX. The overall agreement over the emergent flux in the entire 
wavelength region adopted in our investigation is considered as the evidence for a successful 
implementation of the TLUSTY and the SYNSPEC codes. 

Fig. 3 shows that for a model characterized by T^ff^p = 21000/^ and T^ff^o = 4400 if, 
the SDSS spectrum matches if the distance of the system is about 495 pc and M = 1.5 x 
10^^° y\.Qyr~^. The object J0926+3624 has another component, a spot. Unfortunately, the 
nature and amount of flux contributed by the spot is not known. If the contribution of the 
spot is comparable to that of the disk or the primary, then the present fit of the SDSS data 
at 495 pc implies an upper limit of the flux from the accretion disk. On the other hand 
the observed flux of the primary is a few times less than the model flux of the primary. It 
should be noted that the observed data of the primary is not calibrated for extinction and 
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therefore the actual flux may be slightly higher. Therefore, the distance needed to fit both 
the observed data may be considered as an upper limit. Figure 4 (panel A), however, shows 
that the non-calibrated observed fiux of the primary fits well with a model characterized by 
Tp.ff^p = 18000 -ft" if the di stance is about 465 pc which is consistent with the analysis of 



Copperwheat et al. I (1201 ll ). At a distance d=465 pc, the combined flux also fits well with 
the SDSS data if T^ff^n = 3800 K. If the spot contributes significantly, then a decrease 
in the flux from the accretion disk is needed to fit the SDSS data at this distance. It is 



worth mentioning here that the observation of ICopperwheat et al. I (1201 ll ) doe not indicate 
that the fiux of the spot is comparable to that of primary. The model characterized by 
Tcff,p = 25000 K and T^fj^D = 4000 K is apparently ruled out as the distance needed to fit 
the SDSS data is too high to fit the observed fiux of the primary at the three optical bands. 
A single disk model as described in §2 is considered in all the above comparisons. 

For the range of the donor's effective temperature between 3800K and 4400K and that 
of the primary between 18000K and 25000K, it is found that the flux of the donor dominates 
over that of the primary at wavelength greater than about l.l/i?7i. However, the optimal 
flux of the accretion disk adopted in the present simplifled blackbody model is higher than 
that of both the primary and the donor at all wavelengths. Therefore, future observation 
at infrared wavelengths will enable an estimate of the combined flux of the accretion disk 
and the donor. Since the mass and radius of the donor can be well constrained from the 
evolution model, infrared observation would provide an excellent opportunity to estimate a 
range of the properties of the accretion disk. 

Further, an important flnding in the present investigation is the appearance of emission 
lines in the spectrum of the donor due to the effect of irradiation. Figs. 3 and 4 also show that 
the combined flux does not dilute the emission feature originating from the donor's surface, 
even though the contribution from the disk is higher. This, however is not obvious for any 
combination of the effective temperature of the primary and the donor. As discussed earlier, 
the emission line strength becomes dominant when the temperature difference between the 
donor and the primary increases. In the shorter wavelength region, atomic scattering plays a 
dominant role over absorption. Therefore, the emission lines occur in the near optical region 
only if the irradiation is sufficiently strong so that the scattering albedo is much less than 1. 
However, the temperature-pressure proflles imply strong temperature inversion which indi- 
cates penetration of irradiation energy into the atmosphere even below the region where the 
optical depth is unity. In the region where the optical depth is comparable to unity, photons 
are absorbed and re-emitted in the infrared. Therefore, even with relatively weak irrad i ation , 



emission lines appear in the infrared. On the other hand, JNagel. Ranch fc Werner I (120091 ) 
obtained few emission lines in the spectrum of the disk for an accretion rate higher than 
10~^° M0?/r~^. Therefore, detection of emission lines in the infrared spectrum of the object 



would provide important insight on the atmospheric properties of the donor. 



5. CONCLUSIONS 

Theoretical spectra were calculated for the two white dwarf components in the AM 
CVn system J0926+3624 using the non-LTE atmosphere code TLUSTY and SYNSPEC. A 
range of effective temperatures is adopted based on estimates derived from theoretical binary 
evolutionary considerations for a system formed in the white dwarf channel. Irradiation of 
the cool donor star by the hot primary star is included by assuming complete redistribution 
of energy over the entire atmosphere of the donor star. Observations indicate the presence 
of a spot whose contribution to the total flux is unknown. In the absence of a self-consistent 
non-steady state disk model, a multi-color blackbody spectrum of the disk is considered 
and a maximum disk size is adopted as inferred from observation. Comparison of the total 
flux derived by the present models places upper limits on the distance of the object as well 
as the contribution of the accretion disk to the total flux. The present findings emphasize 
that future infrared observations of the system would provide a good estimation of the 
contribution by the accretion disk to the total flux and hence important insight on the 
thermal properties of the accretion disk. We also predict emission lines in the infrared 
wavelength region, originating from the spectrum of the irradiated donor. The width of 
these lines are dependent on the effective temperature of both the donor and the accretor 
and hence high resolution spectrum may be needed to detect them. If, however, the donor 
is sufficiently cool, as suggested by the evolutionary models, the emission lines could be 
detected easily providing an important diagnostic probe for the donor. 
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quality of this Letter. Thanks are due to C. Deloye for providing the Sloan optical data 
of SDSS J0926+3624 and to F. Allard and D. Homier for providing a set of PHOENIX 
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Fig. 1. — The Temperature-Pressure (T-P) profile of tlie donor star. Tlie numbers near 
tlie lines indicate tlie effective temperature of the non-irradiated donor. The thin solid lines 
represent the T-P profiles of the non-irradiated donor while the thick solid lines represent 
that with irradiation from a primary with T^ff^p = 21000 K and the dashed line with T^ff^p = 
25000 K. The dash-dot lines represent the T-P profiles of the donor with irradiation from a 
primary of T^ff^p = 18000 K. 
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Fig. 2. — The synthetic spectra of the non-irradiated (panel A) and the irradiated (panel B) 
donor star. 7"c//,d is the effective temperature of the donor without irradiation and T^ff.p 
is that of the primary which irradiates the donor. The black lines are the corresponding 
continuum flux. The model with T^ff^D = 4400 K, T^fj^p = 18000 K is also included in panel 
A for a comparison with the non-irradiated case. 
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Fig. 3. — Model fit of the Sloan optical spectrum (red lines) as well as the photometric flux at 
u', g' and r' bands (black error bars). The green line represents the spectrum of the primary 
with effective temperature 21000 K and the black line represents that of the irradiated donor 
with Tcff^D = 4400 K. The dashed line represents the composite blackbody spectrum of the 
accretion disk. The blue line is the combined spectrum of the accretion disk, the primary 
and the donor stars. The results (with different elemental abundances) from the atmosphere 
code PHOENIX are also presented for comparison. 
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Fig. 4. — Same as figure 3 but with different effective temperature of the primary and the 
donor stars. 
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Table 1: The elemental abundances A by number with respect to He adopted for the models. 
Comparison with solar abundances as derived by TLUSTY is presented in column 3. 



Element 


A=N(element)/N(He) 


A/A(Solar) 


He 


1.0 


10.0 


C 


4.79 X 10-5 


0.145 


N 


3.55 X 10-3 


42.7 





2.71 X 10-^ 


0.401 


Ne 


3.83 X 10-^ 


3.19 


Mg 


1.21 X 10-^ 


3.18 


Si 


1.13 X 10-^ 


3.18 


s 


5.16 X 10-5 


2.41 


Fe 


1.01 X 10-^ 


3.2 



